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We propose an efficient linear-scaling time-dependent method for calculating nonlinear 
response function, and study the size effects in non-degenerate two photon absorption 
spectra of Si nanocrystallites by using semi-empirical pseudopotentials. 

1. INTRODUCTION 

The two-photon absorption (TPA) spectrum is an important optical property of semi- 
conductors. It brings valuable information complementary to the one-photon absorption 
spectrum, gives the limits to the optical transparency of materials, and causes laser- 
induced damages. 

Recently, the TPA spectrum of direct-gap nanocrystallites was calculated by Cotter 
et al. QXJ] by using the effective-mass approximation to predict a strong effect of the 
nanocrystallites size. However the experimental result || does not support this result. 
Since then TPA spectra of nanocrystallites has been attracting attention of researchers. 

In this paper, we develop a new algorithm for calculating TPA spectra by using semi- 
empirical local pseudopotentials ||, finite difference method in real space ||], and a linear- 
scaling time- dependent method which has been applied to the calculation of the linear- 
response functions [^|, |7|, §]. This efficient algorithm made it possible, for the first time, 
to calculate the size effect on the TPA spectra of very large nanocrystallites without using 
effective-mass approximation. In the following we show the effectiveness of this algorithm 
by applying it to the TPA spectra of indirect-gap nanocrystallites. Though this result 
does not solve the controversy related to direct-gap nanocrystallites |], 0, it will be an 
important step toward it. 

2. TPA COEFFICIENT 

The non-degenerate TPA coefficient (3i 2 (ui, uo 2 ) describes absorption of the probe light 
with the frequency u\ and the polarization e a in the presence of the excite light with the 
frequency uj 2 and the polarization e^, and, in the transparent region (uji,u 2 < E g ), it is 
related to the third order nonlinear susceptibility Xabba by 

a ( \ 6(2vr) 2 w 1 (3) 

where r] aa (ui) and ^(co^) are the real part of linear refractive index. A simple form of 
^-(3) f or TPA is provided by the second-order time-dependent perturbation theory ||, 
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where the composite matrix element is defined by 

W mf i — ^2 W M — UJi 

The subscripts u and c in the summation run over all valence-band states and all conduction- 
band states, respectively. The subscript m for the intermediate states runs over both 
valence- and conduction-band states. 

It is worth to note that the form Eq. (|2|) looks very similar to that of the linear- 
response function, and that, therefore, we may calculate x by using the linear-scaling 
time-dependent method for the linear-response functions J?]. § . 

x£V)=E — { { ~ ' rVc)(ea ' TCV) - ' rVc){Ga ' TCV) )=( l T dt e+^5B(t)\ (4) 
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where (( • )) indicates the statistical average over random vectors, and SB(t) is the 
response of the system defined by 

8B{t) = 2 lm{$ Ef \e +im {e b ■ r)e~ lHt 9(H - E f ){e a ■ v)\$ Ef ). (5) 

The ket \&E f ) is a random vector projected onto the Fermi occupied states, and 9(H — Ef) 
is a projection operator to extract the Fermi unoccupied states. The imaginary part 
of frequency 7 is introduced to limit the integration time in Eq. (f|) to a finite value 
T = — In 5/7 with 5 being the relative numerical accuracy we need. The main difference 
of from x is the composite matrix elements in place of the dipole matrix elements, 
which we can calculate with the help of the particle source method || ||. 

3. RESULT 
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Figure 1 shows the non-degenerate TPA 
coefficient (3 xx (u\, 0J2) of hydrogenated cu- 
bic Si nanocrystallites of size I — 1 ~ 4 
(nm) as a function of uj\ with a fixed ex- 
cite light frequency, u 2 = 2.4(eV). In the 
calculation, we used the Hamiltonian ma- 
trix discretized into N = L 3 (L = 32 ~ 80) 
cubic meshes in real space, which consists 
of the semi-empirical local pseudopotential 
H the kinetic energy operator in the finite 
difference form . The results were aver- 
aged over 2 — 16 random vectors depending 
on the system size. The energy resolution is 
set to 7 = 200 (m eV), which may not small 
enough to resolve the fine structures in the 
spectra but small enough to study the size 
effects on the magnitude of Imx® 



Fig. 1 Pxx{^i,^2) of Si nanocrystallites 
with u 2 = 2.4(eV). 



The size effects on the TPA coefficient is evident in the figure. The absorption increases 
as the crystallite size increases, and approaches to the bulk value when / 4 (nm). The 
tail extending below the TPA absorption edge 0.8 (eV) in the spectrum is due to the 
Lorentzian distribution with the finite width 7 in the time-dependent calculation. 

4. DISCUSSIONS AND SUMMARY 

The computational cost of our method is O(MN) for the linear response function and the 
non-degenerate TPA coefficient and 0(M 2 N) for the degenerate TPA coefficient. The 
large number M oc T j At oc E max / , ~f <C N is the number of time steps in a time evolu- 
tion or the number of the frequency U\ to be calculated, where E max is the range of the 
eigenenergy. Therefore our method is much more efficient than the 0(N 3 ) diagonaliza- 



tion method |T0[ and the equation of motion methods JTT| whose computational effort is 
0(M 2 N) for the linear -response function and 0(M 3 N) for the third order susceptibility. 

In summary, we have established an efficient linear-scaling time- dependent method for 
nonlinear response function, and studied the size effects on the two photon absorption 
spectra of Si nanocrystallites. Such a large scale calculation has been impossible with 
conventional algorithms. Therefore the present result will be an important step to solve 
the controversy related to direct-gap nanocrystallites JT], . 

The calculation in this article has been done on the supercomputer Fujitsu VPP500 
at RIKEN and NIG. 
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